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(N . Abstract 



In this paper we prove a stable determination of the coefficients of 
the time-harmonic Maxwell equations from local boundary data. The 
argument -due to Isakov- requires some restrictions on the domain. 

Introduction 

Let Q be a bounded Lipschitz domain in the three-dimensional euclidean 

space. Assume the medium, modeled by Q, to be non-homogeneous and 

isotropic. Suppose the electromagnetic properties of Q to be described by 
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the electric permittivity e, the magnetic permeability fj, and the electric con- 
ductivity a. Let E, H denote the electric and magnetic fields, respectively. 



q ■ The time-harmonic Maxwell equations at frequency u > read 

O' 

dH + iue *E = a *E, dE — iuj/j *H = 0, 

whenever the total electric current density is given by a *E. Writing 7 = 
e + ia/u, the time-harmonic Maxwell equations can be expressed as 

dH + iu>y*E = () 

dE-iojfx*H = 0. ^ ' 

It is known that this system may present positive resonant frequencies even 
when the domain is assumed to be of class C 2 (see |23j) or a is assumed to 
vanish (see [15J). This means that, for some positive frequencies, the system 
may have non-trivial solutions of (PQ) with zero boundary conditions. 
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In these notes, we shall study the inverse boundary value problem of 
determining in a stable manner the coefficients /i, e, a by local boundary 
measurements. When setting this problem, the possible existence of resonant 
frequencies makes natural the use of the restricted Cauchy data set as a 
model of non-invasive measurements, instead of using either the restricted 
admittance or impedance maps. Cauchy data sets have been used successfully 
in IU, [21], [22] and ©. 

Let T be a proper non-empty open subset of dfl, the boundary of Q. Let 
v be 1-form defined by v = e(N, .) with N the outward unit vector normal 
to dQ and e the euclidean metric. Given a frequency u > 0, the Cauchy 
data set restricted to T is defined as follows: (T, 5*) 6 C(fi, 7; T) if and only 
if T e TH (T), S E TH(T), and there exists a pair (E,H) E (H(Q; curl)) 2 
solution of ([[]) satisfying *(u A E) = T and *(y A -ff)|r = S. 

In order to quantify the proximity of the restricted Cauchy data sets we 
introduce a pseudo-metric distance which was already used in [6]. 

Definition 1 Let //1, 71 and fi2, 72 be two pairs of coefficients. Consider u> a 
positive frequency and let Cf denote C^Hj^f Y). Let us define the pseudo- 
metric distance between the restricted Cauchy data sets Cp and C 2 as 



S c {Cf, Cj?) = max sup inf || (T J7 Sj) - (T k , S k 



udA &up 1111 \\y±j,jj) — y± k ,o k )\\ TH , T ^ TH , v s. 



\\Tk\\TH (r)—^- 

In order to state our result, we need stable determination of the problem 
on the boundary. Since this has not been proven yet, we shall introduce the 
following definitions. 

Definition 2 Given two constants M, s such that < M, < s < 1/2, 
we shall say that the pair of coefficients /1, 7 is admissible if they satisfy the 
following conditions. 

(i) Uniform ellipticity condition. The coefficients 7, /j, 6 C 1,1 (fi) satisfy 

M" 1 < Re7(x) M" 1 < /x(x); 
for any i6(l. 

(ii) t4 priori bound on the boundary. The following a priori bound holds on 
the boundary 

IMIcMfr) + ll/ i llc°' 1 (r) < M- 



(iii) A priori bound in the interior. The following a priori bounds hold in 
the interior 

Il7llw 2 .°°(n) + IHIw 2 .°°(n) — M, \\i\\h ,2 + s ^i) + llA t llfl-a+*(n) < M- 

Definition 3 Let M, s be the constants given in Definition [2] and let u be 
a positive frequency. We shall say that a pair /i, 7 is in the class of B -stable 
coefficients on V at frequency u if fi, 7 is an admissible pair and there exists 
a modulus of continuity B such that, for any other admissible pair ji, 7, one 
has 

117 ~ Tllc°.i(r) + IIa* - Allco.i(T) ^ 5 f fc(C r , Cr) j , 
l|V(7-7)ll LOO( r ;C 3) + 11^(^-/2)11^(^3) < B (6 C (C T , Cr-)) . 

Here Cr, Cr are the Cauchy data sets associated to the pairs fi, 7 and /2, 7, 
respectively. 

The first idea in our argument is to construct special solutions vanishing 
on dQ \ T, the inaccessible part of the boundary. In [10] Isakov proposed a 
reflection argument which allows to construct solutions for the conductivity 
equation with the desired behavior on the boundary. This argument was 
extended in [5] to the time-harmonic Maxwell equation. 

In order to carry out Isakov's approach it seems to be necessary to assume 
some geometrical restrictions about the domain, namely, the inaccessible 
part is supposed to be either part of a plane or part of a sphere. Despite 
this restriction, the method allows to prove the following result from local 
boundary data. 

Theorem 1 Let U be either a suitable partially flat domain or a suitable 
partially spherical domain and let wbea positive frequency. Consider ^1,71 
and fj,2, 72 any two pairs in the class of S-stable coefficients on T at frequency 
cu, with B satisfying \r\ < B(\r\) for all \r\ < 1. Assume that d^Hj = d^j = 
on dll \ T with j = 1,2. Then, there exists a constant C = C(M) such 
that the following estimate holds 

||7i - 72||jri(t7) + ll^i ~ ^\\ H i {u) < C [log B(6c(C T ,C^))\' x , 

for some constant A such that < A < s 2 /3. Here C r , C^ are the restricted 
Cauchy data sets associated to the pairs yUi,7i and yU.2,72, respectively. 

The exact meaning of suitable is explained in Subsection 11.41 



As in the inverse conductivity problem, it should be possible to prove 
that any admissible pair is in the class of Holder-stable coefficients on T for 
any frequency u>, that is, with -B(|r|) = \r\ a for < a < 1. Note that we 
have obtained the same kind of stability as in the global data case (see [6]). 

From the point of view of applications it might be useful to suppose 
the coefficients to be equal on the accessible part of the boundary. In this 
particular case we get the following corollary. 

Corollary 1 Let U be either a suitable partially flat domain or a suitable 
partially spherical domain and let wbea positive frequency. Consider /ii,7i 
and //2,72 any two pairs of admissible coefficients. Assume that 

Mi|r = A^lr, d&n^v = d&^lr, 7i|r = 72k, <9^7i|r = d x j^ 2 \r 

and djsffik = QnIu — on dU \ T with j = 1, 2, 3 and k — 1, 2. Then, there 
exists a constant C = C(M) such that the following estimate holds 

II71 -l2\\ H i {u) + Whi-IMiWhiqu) < C|log5c(Cr,Cr)r A , 
for some constant A such that < A < s 2 /3. 

Furthermore, if we follow the proof of Theorem Q] one can state the fol- 
lowing uniqueness result. 

Theorem 2 Let U be either a suitable partially flat domain or a suitable 
partially spherical domain and let wbea positive frequency. Consider /ii,7i 
and /i2,72 in C 1,l (U) such that 

Mi|r = ^Ir, d x j^i\ T = d xJ fi 2 \r, 7i|r = 72k, <9^7i|r = d xi j 2 \r, 

with j = 1, 2, 3. If additionally Cf = C£ and &N^k = d^lk = on dU \ T 
with k — 1, 2, then 

A*i = A*2» 7i = 72 

inU. 

As in the inverse conductivity problem, it should be possible to prove that 
the coefficients are equal on the accessible part of the boundary T whenever 

Op — Lyp. 

The problem of determining the electromagnetic coefficients by data taken 
on the entire boundary has been studied by several authors. The unique 
recovery of C 3 -coefficients 7 and \x from boundary data was proved in [18], 
and later simplified in [20]. Boundary determination results were given in [12] 



in the case that the boundary is smooth. The more general chiral media was 
studied in [TB]. For a slightly more general approach and more background 
information, see also the review article |19| . 

The inverse problem of determining the electromagnetic coefficients from 
partial data has been much less considered. As far as the author knows the 
only work in that direction is [5]. 

Two different approaches have been used to attack the inverse conduc- 
tivity problem from partial boundary data. The first one was proposed in 
[3] and generalized in [14] . In [7j, this method was used to give a log- log- 
stable determination in the framework of [4J. In this approach there are not 
any strong geometrical restriction about the domain but the partial measure- 
ments have to be taken in the whole boundary. Getting an optimal stability 
(i. e. a stability with a log-type modulus of continuity) in the context of [4J 
may be difficult and the stability for [13] is an open question. The second 
approach for partial data was proposed in [10] and the optimal stable deter- 
mination was stated in [8]. As we have already mentioned, this argument 
requires a strong restriction on the domain. However, the measurements are 
localized on the accessible part of the boundary and it is possible to get the 
optimal stable determination. These two facts are very important from the 
point of view of applications. For instance, Alessandrini and Vessella proved 
in [2] that a logarithmic estimate yields Lipschitz stability for some finite 
dimensional spaces of conductivities. 

These two approaches have been extended to systems. In [22] Salo and 
Tzou followed the spirit of [H] to prove uniqueness in the context of Dirac's 
equation. Isakov's argument was extended in [5] to Maxwell's equations. 
The proof given here takes some ideas from [13] and it turns out to be more 
convenient and useful for us than the proof given in [5]. In fact, it avoids 
the long computations made there to prove the thesis of Theorem [2] and 
it allows to relax the hypothesis about the domain and the smoothness of 
the coefficients. In [5] the domain was assumed to be of class C 1 ' 1 and the 
coefficients were assumed to be C 4 . Besides, a technical hypothesis about 
the extension of the coefficients had to be supposed. 

An overview of the paper is the following. It has three sections. In the 
first one, we give some preliminaries about the functional spaces used. In the 
second section we prove our results when U is partially flat. To achieve this, 
we use a reflection argument to construct special solutions vanishing on the 
inaccessible part of the boundary. In the third section we connect the flat 
case with the spherical one by means of the Kelvin transform. 
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1 Preliminaries 

Let E be the three-dimensional euclidean point space and let its tangent 
bundle be denoted by TE. Let TE be the module of smooth vector fields 
over the real smooth functions C°°(E; R) and define 

XE = {u + iv:u,ve TE}. 

The elements of XE will be called complex vector fields. Let the bundle of 
alternating tensors be denoted by A k TE with k — 0, 1, 2, 3. Let A k E be the 
vector space of differential /c-forms and define 

A k E = {oo + ir] : oo,7] eA k E}. 

The elements of A fc E will be called complex fc-forms. Recall that 0-forms are 
smooth functions by definition. As it is usual, d and A denote the exterior 
derivative operator and the exterior product of forms, respectively. 

The euclidean metric e induces a volume element denoted by dV, a dis- 
tance denoted by d e and a point- wise inner product denoted by (u, if) for any 
uu,r] G A k E with k = 0,1,2,3. Recall that the Hodge star operator is the 
unique bundle map * : A k TE — y A 3-A TE satisfying 

w A *rj — (uu, i]) dV. 

Moreover, * *u = u. Let us define \r]\ 2 = (77,77). 

The formal adjoint of d will be denoted by 5 and it can be expressed by 

5r] = (— l) k *d*i] 

for 77 G A k . Let us define the laplacian on /c-forms as —A := Sd + dS. 

We also recall that we can identify vectors and 1-forms by means of the 
metric, that is, 

ueTEi — y 77 = e(u, .) e A X E. (2) 

If u G TE, its corresponding 1-form will be denoted by u . However if the 
difference is clear by the context it will be denoted by u. On the other hand, 



if v G ^4 1 E, its corresponding vector field will be denoted by vK As before, 
this notation will be used whenever the context is not clear. 

Finally, for any / G C°°(E; R) and any u, v G 7TE, u ■ v = e(u, v) denotes 
the point-wise inner product, uxv — *(w b Av b )^ denotes the point-wise cross 
product and V/ = (d/)", V • u = — Su and Vxm = (*G?w b ) tl stand for the 
gradient, divergence and curl, respectively. 

1.1 The functional spaces 

Along these notes we shall say that a domain Q is Lipschitz if its boundary 
dQ is locally the graph of a Lipschitz function. Additionally, N denotes the 
outward unit vector normal to dCl and v := iV b is its corresponding 1-form. 

In order to perform the proofs of the results stated in the introduction we 
require some standard functional spaces: H S (K),H S (Q),H^(Q) with s G R 
denote the potential Sobolev spaces based in L 2 ; W 2 '°°(Q) stands for the 
Sobolev space with two derivaties in L°°; B s (dQ) with < \s\ < 1 denotes 
the Besov spaces B* (d£V) with p = q = 2. A quite complete description of 
these spaces can be found in [IT] . 

Additionally, when working with Maxwell's equations other non-standard 
Sobolev and Besov spaces turn to be useful. Those are mainly, H(Q;div), 
if (Q; curl) and TH(dQ). The first one corresponds to the fields or 1-forms 
in L 2 with divergence in L 2 . In this space, it makes sense the normal traces 
as elements of B~ l l 2 (dVL). The space H (£); curl) corresponds to the fields 
or 1-forms in L 2 whose curl is in L 2 . The tangential traces of elements of 
H(n;cur\) make sense as elements of B-^dtyC 3 ) or B- x ^ 2 (dSl;A 1 TE). 
The space TH(dQ) is defined as the space of tangential traces of H(Q; curl) 
and we have that 

HHIs-i/a^OiC 3 ) + ll Divu7 lls- 1 /2(an) ^ Ci IMIta^o) ' ( 3 ) 

WMIthou) < °2 (l|w|| B -i/ a(8 n ; c3) + W Divw \\B-^(dn)) > ( 4 ) 
if w is a vector field, or 

ls- 1 /2(an;A 1 TE) + lPivw|| B _i/2( 9n ) < Ci ||u>|| Tjf/(an) , (5) 



w 



ll W llra(ai) ^ ^2 ^||w||B-Va(«l;AlTE) + \\ B ' W W \\ B-^(dQ)) > ( 6 ) 

if w is a 1-form. Here Div stands for the surface divergence. Recall that the 
surface divergence of w G TH(dQ) makes sense, and it can be defined as an 
element of B l l 2 (dVL) in the following way: 

Divu> = — JV'(Vxu), Divu; = — (u, *du) 



where u G H(Q; curl) and Nxu = w, if u is a vector field, or *(u Au) = w, 
if m is a 1-form. 

Finally, we recall some key points. For any u G H (Q; div) and any g G 
B l ' 2 {dVt) we have 

(JV-%)= f(V-u)JdV+ [ u-VfdV, (7) 

where w is a vector field, / G if 1 (fi) and /|an = g. We also have 

((i/,u)|<7) = - f(6u)JdV+ f(u,df)dV, (8) 



where w is a 1-form, / G if 1 (f2) and /|gn = g. 
The maps 

iVx. : TH(dVL) — ► (TH(dVL))\ *{v A .) : TH(dVt) — ► (TH(dVt)Y 

are isomorphisms. In particular, 



(Vxu)-Ud^= / u-(Vxv)dV-(JVx4/Vx(iVxv)), (9) 



for vector fields, and 



*du,v)dV = (u,*dv)dV - (*(v Au)\*(v A*(v Av))} , (10) 

in 

for 1-forms. 

A detailed exposition of the collected facts can be seen in [T7] and [B]. 

1.2 Some remarks on the boundary 

Along this section, fi denotes any bounded Lipschitz domain, dQ denotes its 
boundary and T stands for a proper non-empty open subset of dQ. Moreover, 
.|r and .|p denote the restriction to T and T, respectively. 

Definition 4 For < s < 1, define the space B s (T) as 

B S (T) = {/| r : / G B s (^)}, 

with the norm 

I|0|Ib-(D = inf {ll/ll^(aQ) : /Ir = 9}- 



On the other hand, for < \s\ < 1, define 

B s (T) = {feB s (dQ):snppfcT}, 

with norm 

11/ II.Bg(r) = ll/lls s (an) • 
Finally, for < s < 1, define the space B~ S (T) as the dual of Bq(T), that is, 

B- S (T) = (B S (T)Y. 

Note that, for < s < 1, Bq S (T) is the dual space of B s (T), that is, 

B S (T) = (B S (T)Y. 

Lemma 5 Let s, e be such that < s < 1 and < e < 1 — s. Then there 
exists a constant C(s, e) > such that, 

(a) for any g G C°' s+e (dtt) and any / G B s (dty, 

\\9f\\B s (dn) — C lbllc°. s + e (9a) H/llB s (on) > I 11 ) 

(b) for any g : dVt — > C with g G C°' s+e (T) and any / G B S (T), 

||<7/||bj(T) ^ C IMIc°.H"(T) H/Hbs(T) 5 ( 12 ) 

(c) for any # : <9fi — ► C with # G C°> s+t (T) and any / G B s (F), 

\\9\rf\\Bs(r) — C \\9\\co, s +e(p) \\f\\B*(r) ■ (^) 

Remark: The constant C(s, e) given here blows up when e becomes 
small. 

Proof: (a) was proven in [6]. (b) and (c) follow easily from (a) taking 
g G C°' s+e (dQ) an extension of g\f such that 

ll^llc°> s + e (dn) — ^ IMIc'^+^r) • 
Regarding to extensions from a closed subset of E, see [21] . □ 

Definition 6 Define the space TH(T) as 

TH(T) = {w\ r : w G TH(dn)}, 



with the norm 

\\4th{t) = ' mi {\H\ T H( 9 n) : w \r = A- 
On the other hand, define 

TH (T) = {we TH(dn) : suppw C f}, 

with norm 

IMIthoIT) = IMIth(<9£7) • 

Lemma 7 Let N be the outward unit vector normal to dQ and let v be its 
associated 1-form. Then 

NxTH (T) = (TH(T))*, *{y A TH (T)) = (TH(T))*. 

Proof: Here we prove the first identity. The second one follows by the 
correspondence between vector fields and 1-forms. 

Let I : TH(T) — > C be a bounded linear functional, we can construct 
another functional I : TH{dVt) — ¥ C defined by l(w) = l(w\r), for any w G 
TH(dQ). Since / is linear, bounded and 



(TH(dn)r - l|Z| W))* 



there exists z G TH(dQ) such that (Nxz\w) = l(w) = l(w\r) with 

lpllTif(af2) — lrll(Ti?(r))* ■ 
It is clear that suppiVx,2 C T, hence z G TH (T) and 

IMIt_Ho(T) — lrll(Ti?(r))* • 

Conversely, given z G TH (T) we can define / : TH(T) — > C as l(w) = 
(Nxz\w), for any w G TH(T) and w G TH(dQ) such that w\t = w. It is 
well-defined since suppiVxz C T. Moreover, 

\l(w)\ < \\z\\ THo ^ \\w\\ TH ( dn ^ , 

which implies 

\l{w)\ < \\z\\ THo ^ \\w\\th(t)- 
Therefore, / is a bounded linear operator with norm 

lnl(TH(r))* — IMl77f (r) ■ 

D 

10 



Lemma 8 There exists a positive constant C such that: 

(a) For any w G TH(dQ) and any / G C 0,1 (dQ), one has that 

Wf w \\TH(an) ^ C ll/llc - 1 ^) H w llTH(9n) • ( 14 ) 

(b) For any iu G TH (dtt) and any / : «9fi — ► C such that / G C 0>1 (r), 
one has that 

\\fw\\ THo(r) < C ||/|| c o,i ( r) \H\ T H (r) ■ ( 15 ) 

Proof: (a) was proven in [5J. (b) follows easily from (a) taking and 
extension / of /|p such that / G C 0,1 (dQ) and satisfying 

Again, regarding to extensions from a closed subset of E, see |24j . □ 

Lemma 9 The following items hold: 

(a) If w E TH (T), then Divw G 5 ~ 1/2 (r) and 

||Div™|| B -i /2(r) < C \\w\\ THo{v) ; (16) 

(b) If z G TH(ffl), then, for any / G Sj /2 (r) and 5 G TH(dn) such that 
z|r = ^|r, one has ((Divz)| r |/) = (Divz|/) and 

||(Divz)| r || B -i/ 3(r) < C \\z\ T \\ TH{v) ■ (17) 

Proof: It is easy to check both items. 

(a) Divw is well-defined and belongs to B~ l l 2 (dVt). It remains to prove 
that supp Divw C V. In order to verify this last point, we just need to 
have in mind the following facts: if / G H 1 ^), then NxVf G TH(d£l); 
moreover suppiVxV/ C supp/|ao. Indeed, supp Divw C T since 

(Dww\f\ dn ) = (w\Nx(NxVf)) . 

The estimate is now immediate using either (EJ) or (J5J). 

(b) By an analogous argument to the one given in (a), we have that if 
z\t — z \r then (Divz)| r = (Div^)|r- Hence the identity follows. The 
estimate is a consequence of the identity and (J3]) or (J^J). 

□ 
11 



1.3 Maxwell's system as a Schodinger equation 



In this section we shall transform Maxwell's equations into a Schrodinger- 
type equation. The idea of this transformation was already introduced in 

There is a well known process which allows us to transform Maxwell's 
equations into a Schrodinger-type equation. In order to do so we require some 
extra smoothness of the coefficients, namely //, 7 € C l,1 (Q). The first step in 
this process is to augment the Maxwell system with two scalar equations: 



0. 



8(jE) = 0, 8{nH) 

The information coded in these scalar equations was already present in the 
initial system. In order to check this, it is enough to take *d in each equation 
in flU). 

Next, we introduce a new system inspired in the four mentioned equations. 
This new system reads as 

5(jE) + iwyfih = 

— 7~ x <i(7e) + *dE — iu>/j,H = , . 

H^dfah) + *dH + iwyE = 0. 

The new terms preserve the physical units of measure of the original four 
equations. Choosing euclidean coordinates, the new system -called hence- 
forth augmented system- can be written in vector field notation as it follows 



/ 


D 


D- \ 

-Dx 


D- 

\ D Dx 




J 



( Ufl 



Ufll 3 



DP- 



\DP 



Da- \ 



Da 



wy 



W7/3 / 



f h\ 
H_ 
e 



0. 



where a = log 7, (3 = log/i, Ij is (J x j)-identity matrix, with j e N and 

D- = - ( d^i d x 2 d x 3 ) , 



- ( d x i 
1 v 



D 




i 



d x 3 
-d x 2 



-d x 3 



Or 



d x 2 



In a much more compact manner we shall express the augmented system as 
(P + V)X = 0, where 



P 



( 


D 


D- \ 

-Dx 


D- 
\D Dx 




) 



( un 



V 



u/jlI 3 



Dj3- 



\D/3 



Da- \ 



Da 



CU7 



W7/3 / 
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Note that E, H is a solution for Maxwell's equations, if and only if, X 1 = 
( h H * I e E t J is a solution for the augmented system and the scalar fields 
e, /i vanish. 

The next step is to rescale the augmented system, that is 



(P + V) ( 


' ^ 1/2 h 






V 


T 1/2 h 


where 


W = Kh + \ 


( 







Y 



1 



■wl 



// 



W 



(P + W)Y, 



Da- \ 

Da Dax 



D/3- 
\ Dp -D/3x J 

with k = bj^l 2 ^ 1 ! 2 . We shall call 

(P + W)Y = 

the reseated system. 

The advantage of rescaling is that 

= (P + W)(P - W l )Z = (-A/ 8 + Q)Z, 
= (P - W l ){P + W)Z' = (-AJ 8 + q')z', 
= (P + VT)(P - W)Z = (-A/ 8 + Q)Z, 



(19) 



(20) 
(21) 

(22) 



where Q,Q',Q are zeroth-order terms. Here W f denotes the transposed of 
W and W* stands for W l . No first order terms appear in (J2DD, (ED) and (j2"2"|) . 
giving as a result a Schrodinger-type equation. Mind 



Q = -PW 1 + WP- WW*. 



(23) 



Note that if Z is a solution for ([20]) in fi, then Y = (P - W*)Z is a solution 
for the rescaled system in $7, hence 



X 



/' 



-V2/ 4 



7" 



^ 



y 



is a solution for the augmented system. In the same manner, if Z is a solution 
for Q22D, then Y = (P - PF)^ is a solution for (P + W*)y = in fi. 
For later uses, 



/ A 



^ 2 



a 



2V 2 a - AaP 



V 



\ 



A(3 



+ 



2V 2 /3- A/3/3 / 
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«* 



\{Da-Da))I A 



2Dk- 



\ 2Dk 

( A/3 



2Dk- \ 



2Dk 



«* 



\{Dp-DP))h 



(24) 



/ 



Q / = -i 



2V 2 /3 - A/3/3 



V 



\ 



Aq 



K 2 + MD(5-Df3))h 



\ 



-2Dkx 



2V 2 a-AaI 3 j 
2Dkx 



[K 2 + \{Da-Da))h 



I 



(25) 



and 



Q 



1 


/ 


-A/3 


-2V 2 /3 + A/3/3 




\ 


2 


V 






-Aa 


-2V 2 a + Aa/ 3 j 



( <-± I 



K 2 -\(D0.D0))I A 



\ 



2Dkx 



\ 



-2Dkx 



k 2 - \(Da-Da))h 



(26) 



/ 



with V 2 / = (% M /)J >fc=1 . 

In order to make as concise as possible the presentation of our proofs, we 
introduce some additional notation. Let Y, Z be in the form 



Y=(f 1 (u 1 )'!/ 2 (« 2 )*)*, Z=((/ 1 W|j 2 (v 2 Y)\ 



define 



(Y\Z) = V ( f f j gi dV + f u j -vi dV) , 

=1 V^n in / 

(y\Z) d n = J2([ f 3 &dA+ [ 
■ =1 V^an ./a 



u 3 -vidA . 



In the first identity we are assuming / J ', <7 5 ' G C°°(f2) and w- 7 , f- 7 G ^E|jy with 
j = 1,2, while in the second identity p, gi G C°°(dQ) and vP ,vi G A'Ejgn 
with j = 1,2. The following integration by parts holds 



(py|z) = (p n y\z 



dn) dn 



(Y\PZ). 



14 



Here, when A is a (possibly complex) vector field we denote 



/ 


A 


A- \ 

-Ax 


A- 

\ A Ax 




) 



(27) 



Finally, for elements Y in the form given above we define, for \s\ > 0, 

livii — V^ i II pII -i- \\v j \\ i 

W 1 \\H°(n;y) - Z^ \\\ J ll-ff s (n) ^ II" \\h°(ci;C s )J > 



i=i,2 



and 



mix» ( noo = E(ll/ i ii— + 11^' 



Il 2 (h) 



L 2 (Q;C 3 ) / ' 



i=i,2 



On the other hand, we define, for < \s\ < 1, 



IVII — V^ I II Pll -l- 11?/-? I 



i=i,2 



is s (an) 



_B s (df2;C 3 ) 



and 



ivil - V^ i II f?ll -i- IU/ J 'I 

I 1 llL2(an ; y) - 2^ III-' H>-2/»oA-t I" 



iL 2 (an) 



L 2 (d£7;C 3 ) / ' 



3=1,2 



1.4 About the geometry of U 

In order to make precise the geometrical restrictions assumed in Theorem [TJ 
Corollary [1] and Theorem El we give the following definitions. 

Definition 10 We shall say that a bounded Lipschitz domain U C E is 
partially flat if there exists a plane g C E and some euclidean coordinates £ 
such that, 

(i) g = { x e M. 3 : x 3 = 0}, 

(ii) f/c{iei 3 :i 3 < 0}, 

(iii) T := mt q (dU n g) ^ 0. 

We shall say that a bounded Lipschitz domain U C E is partially spherical 
if there exist a point Qo £ E, r > and some euclidean coordinates £ such 
that 

(i) Q = y and U C .6(2/0; r ) := {y eR 3 :\y- y \ < r }, 
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(ii) T := mt S ( yo ;ro)(dU n S(y ; r )) 7^ wh ere S(y ; r ) := <9£(?/o; r ), 
(iii) e S(y ; r ) but i U. 
In the two previous cases, we denote V := dU \ F . 

Definition 11 We shall say that a partially flat domain U is suitable if its 
symmetric extension with respect to q -that is Q := U U To U TZ(U)- is also 
Lipschitz. Here 1Z denotes the reflection with respect to q and it is defined 
as (x^x^x 3 ) 1 — > (x 1 ,x 2 ,—x 3 ). 

In addition, we shall say that a partially spherical domain U is suitable 
if its inversion with respect to 5(0; 2r ) -that is Q := ]C(U)- is a suitable 
partially flat domain. Here K denotes the inversion with respect to S(0; 2ro) 
and it is defined as y 1 — > r\/\y\ 2 y with r\ = 2r . 

We have to restrict ourselves to these suitable domains because we need 
to make an extension of the coefficients preserving their smoothness (see 
Subsection 12.31) . 

2 The domain U is partially flat 

Along this section we assume U to be a suitable partially flat domain and 
we follow the notation in Definition [10] and Definition [TTJ 

2.1 Maxwell's system and the reflection map 

Let the coefficients //, 7 be such that //, 7 G C l,1 (U) with x 3h\y = d x 3j\r = 
and set ft, 7 : Q — > C two smooth extensions of /i and 7 defined as 

/^(x 1 ,^ 2 ,^ 3 ) = ^(x 1 , x 2 , — |x 3 |), 7(x x , x 2 ,x 3 ) = 7(x 1 , x 2 , — |x 3 |), 

for any x & Q. Note that the hypothesis d x 3fi\ ro = <9 x 37| ro = allows us to 
keep the smoothness when extending. 
Consider the system 



V x H + iuo^E = 
VxE - iujpiH = 

in Q. The push-forward of the reflection map TZ reads 

1 
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(2* 



Let / be a smooth function and u, v! two vector fields on E. Let g,v,v' 
denote the function and the vector fields given by 

g(x) := f(R(x)), v x := U^u n{x) , v' x := TZ^u n{x) . 

It is a straight forward computation to check that 

(V-v)(x) = (V-u)(n(x)), (29) 

(Vg) x = K.(Vf) nx) , (30) 

(Vxv) x = -K*(Vxu) n{x) , (31) 

(v x v ') x = -11* (u x «') K (x) • (32) 

On the other hand, let a be a smooth function defined in {x G M 3 : x 3 < 0} 

and set a, the extension of a to E, defined as a(x 1 1 x 2 1 x' i ) = a(x 1 ,x 2 , — \x 3 \). 
Then 

(Va) x = 7^(Va) w(a0 . (33) 

Lemma 12 Given 

Y= ( ji^H^O f^E 1 )\ 

such that E,H G i/(fi;curl) is a solution of (128 p in Q, one has that £? — 
E } H -H } with 

-E-z := lZ*En( x ), H x := —1Z*Hn( x y, 

is also a solution of (|28l) in Q satisfying 

Nx(E-E)\ ro = 0. (34) 



Proof: Let E, H be a solution of (|28l) in fi. It is an immediate conse- 
quence of (13T|) and the definition of /i,7 in O that E, H is also a solution 
for (1281) in Q. Further, from the weak definition of tangential trace one can 
derive that_iV x (E - E)\ Fo = 0. Indeed, let w G B 1 / 2 (dU; C 3 ) such that 
supp w C r and consider v G H 1 ^; C 3 ) such that v\qu = w, then 



NxE-NxE 



\ = J (VxE-VxE)-vdV - J (E-E)-VxvdV 

I VxE-v- E-VxvdV + I {V x E) n -V~v + En-Tl^xvdV 
Ju Ju 



WxE-v-E-WxvdV + / (Vx£)-^w K -£-(Vx^ R )(il/ 
c/ Jn{u) 



WxE-u- E-WxudV = (NxE\u) dn = 0. 



Here we have used fl3Tj) twice, and the fact that w, defined as v in [/ and as 
K*v n in ft(E7), belongs to H 1 ^; C 3 ) and u| an = 0. □ 
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Lemma 13 Given 

Y={f l u'\p M 2 ) 

solution of (P + W*)Y = in Q with p e H\Q) and u j G H(Q; curl), one 
has that Y - Y is also a solution of (P + W*)(Y - Y) = in H Here W 
denotes the matrix (TH?]) for coefficients /2, 7 and Y x := JY-j^m with 



/l 



Y x = ( /i(a;) «i I f\x) u\ ) 
Additionally, 



J :- 



-11* 



\ 



\ 



n* J 



(f - h\r = 0, iVx(n 2 -n 2 )|r 
Proof: The first part of the lemma follows from 

(PY) X = J(PY) n{x) 



and 



(35) 

(36) 

(W*Y) X = J{W*Y) n{x) . (37) 

The identity fl36|) is a consequence of f )29|) . fl30|) and fl3T| . The identity f )37|) 
follows from (|33|) and fl32j). 

Additionally, (J 1 - j 1 )]^ = since f 1 G H 1 ^) and iVx (m 2 - u 2 )\ ro = 
as we showed in the proof of Lemma [T2J □ 

2.2 Relating the boundary measurements with the co- 
efficients in the interior 

Lemma 14 Let ^-,7, belong to C 0,1 (U). Then, for any Y\ given as in the 
hypothesis of Lemma [T21 with coefficients /ii,7i and any Y 2 given as in the 
hypothesis of Lemma [TBI with coefficients /a2,72, one has that the elements 
£1 = Ei — Pi, H\ = Hi — Hi and y 2 = Y 2 — Y 2 , expressed in the form 

*i = (o ti /2 m\o 7 1 1/2 ^) i ^ 2 = ( /1 {ul)t \ f2 W> 

satisfy the following estimate 



{Yi\P72) n -{PYi\72) n \< 



<C8 c {ClC*) 



-1/2 
M 2 



lhlr|| Bl/2(r) + 



C ' 1 ^) 



72 



1/2 



— pi r TH/n ~r 
c*o,i(r) A n \ L ) 
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-1/2 
72 



c°,i(r) " UB o ( l > 



1/2 



- \\ z 2 
CM(r) 



\th (t) j ll^ x ' E illr J Ffo(r) 



+C 



-1/2 -1/2 

Mi - lh 



cm(d 



|<?2|r|| B i/2 (r) 



+ 



1/2 1/2 

7i -72 



co.ifr) 



Mr|li\ff(r) + 



-1/2 -1/2 

7i -7 2 



+ 



1/2 1/2 

Hi -f*2 



c* - 1 ^) 



\Z>2 



lrflb(T)J (ll^x^il 



c o,i(r) " IIB o ( r ) 



TH (T) 



+ \\NxXi\r\ 



TH(T) 



Here gx , 32 E B l l 2 (dU) stand for £1 = f l \ du , 32 = f 2 \ du and z U z 2 E TH(dU) 
stand for z\ = Nxu 1 , z 2 = Nxu 2 . Here C 3 T = C([ij, jj\ T) with j = 1, 2. 



Proof: It is easy to check, using (1361) that 

where 9^ = Y"i — Yi. Let x be 

x = ( /if^|0 7 2 /a E|)', 
with £2, #2 £ H(U; curl) an arbitrary solution of 

V X ^"2 + i^l2^2 =0, V x e 2 — iOJ^2^2 = 



(3? 



in [/ and satisfying supp iVx£ 2 C 17 Since (P+W^)^ = and (P+W 2 )l = 
in U, one has that (xlP^)^ = {P^\f2)xji hence 

WPotiu - {PirMv = (* - xIP^ - (p(^ - ^)|r 2 ) c/ • 

On the other hand, we have, using (J2J) and (jHJ) , that 

(fTi - x|Pr 2 )[; - (P(fXi - l)\j 2 ) u = 
= ilN-^Hi - 112X2) H2 1 32) +i(N'(fi 1 y{i) (>7 1/2 - fi 2 1/2 ) 32 



+z(iV-(7i£i -7 2 £ 2 y 



-1/2 
7 2 3\ 



-ilNxirtx-Hi) 



}+^(iV-( 7l x 1 )|(7r 1/2 -72 _1/ V) 
)|iVx( / if, 2 ))-z(iVx^ 1 |iVx((^ /2 - / if), 2 )} 

-t (iVx ( 72 1/2 iVx ( El - £ 2 ))|*i) - % (Nx (( 7l 1/2 - 72 1/2 )iVx £l )|^) . 

Furthermore, from the Maxwell's equations one deduces that 



N-dj'Ej) = — Div(iVx^), iV-(^) 



/lJ 



lijj 



Div (Nx<Ej) 



19 



for j = 1,2. Hence, 



GXil-P^r/ - (-Prilf^X 



w 



-(Div(iVX£i -iVx£ 2 ) 



l''2 1/2 S2) (DiviVxEi 



/ -1/2 

Ol - ^2 



l^ 



;^2 



+^ (Div^x^ - iVx^ 2 )| 7 - 1/2 ^) + I ^DiviVx^ 1 |(7 1 - 1/2 - 7 - 1/2 )^) 

Nx{nl /2 z 2 )} - i (nxX^Nx ((fi{ /2 - fil /2 )z 2 ) 



-i(Nx!Hi- Nxti 2 



-^(iVx( 72 1/2 (iVx£ 1 -iVx£ 2 ))|^ 1 )-^iVx(( 7l 1/2 -7 2 1/2 )iVx£ 1 )|^ 1 ). 

Let us denote Nxt.j = Tj, Nx9(j\ T = Sj, then by using the appropriate 
dualities, the boundary conditions (IMj) . ( Bo'j) and the estimates (!T2"j) . (IT5|) . 
(USD, (USD and flUD we get 



I (Vil^aJn - ( PY M 
< C [\\Ti - T 2 \\ THo(r) + || 5i - S 2 \\ TH{r) 

\ z ^\\th(t) + 



< 



1*2 



1/2 



+ 



72 



1/2 



c^md 



-1/2 
72 



1/2 



co.1^ W^ocn) + c ( 



-1/2 



-1/2 
/"2 



c o,i ( r)^ 2|r|lsl/2 ( r ) 
WIb 1 / 2 (t) 

lklr|| Bl/2(r) 



C.^r) 



+ 



1/2 1/2 

7i -7 2 



c . 1 ^) 



\ z Mn\TH{V) 



-1/2 -1/2 

7i -7 2 



- \\3i R^/n 
co,i( r ) ii ii-Kq I. 1 J 



+ 



1/2 1/2 



C - 1 ^) 



l z 2|lr^o(r) 



l^il 



T^o(r) 



+ l|Si| 



TH(r) 



This estimate holds for all (T 2 , S 2 ) G Cp, since £2,^2 was chosen to be an 
arbitrary solution of (I38I) in £/ satisfying supp iVx£ 2 c T. Finally, the wanted 
estimate is a consequence of Definition [TJ □ 

Proposition 15 Let 71, [i\ and 72, /x 2 be in the class of S-stable coefficients 
on T at frequency u>. Then, there exists a constant C(M) such that, for any 
Z\ E H 1 ^; y) satisfying Y\ = (P — W{)Z\ with Y\ as in Lemma [Til and any 
j 2 G H l (Q; y) as in Lemma dU one has 



((Qi-Q2)Z 1 \or 2 ) a \ <CB(5 c (C l Tl Cl)) ||ZiL W) INI^ 



+ CB(5c(C l T ,C 2 T ))(\\>E l 



Ijf(l7;curl) + ll-^lliir^curl) 



Hnjy) 



v 1 II/ 1 !! -i- II Ml -i- ll/ 2 ll -i- II 2 II 

\lK llfl-i(I/) "^ ll U llff(t/;curl) ~^ II* llif 1 ^) ~^ II" \\H(V\cox\)} ' 

Here Qj is the matrix (T23J) associates to /ij, 7,- with j = 1,2. 
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Proof: From 



one has 



((Qi - Q 2 )^ 2 ) Q = - (P{W{ - w^)z 1 \^ 2 ) n 

+ ((W x - W 2 )PZ 1 \y 2 ) n - ((W X W{ - W 2 W t 2 )Z 1 \y 2 ) n 
= ((W{ - W 2 t )Z 1 \P N y 2 ) m - (W t 1 Z 1 \Py 2 ) n + (w 2 t z 1 \p? 2 ) n 

+ (w 1 (p - wi)z 1 \or 2 ) Q - (pz 1 \w*? 2 ) n + (w 2 %|w 2 v 2 ) n 

= (iWl - Wl)Z x \P N y 2 ) m + ((P - W t 1 )Z 1 \P^ 2 ) n + (W X (P - Wl)Z x \<y 2 ) u 
= ((W[ - W 2 t )Z 1 \P N ^ 2 ) m + (Y 1 \Pf 2 ) il - (PY x \K) n . 

In order to get the penultimate identity, we used twice that (P + W 2 *)y 2 = 0, 
while to get the last one, we used that Y x = (P—W{)Z X and that (P-\-W\)Y\ = 
0. 

It is a straight forward computation to check the next estimate 

| ((W[ - W 2 t )Z 1 \P N Y 2 ) m | < C (\\ Kl - n 2 \\ Loc{r) + \\V(fa - /3 2 )|| Loo(r;C a ) 

+ ||V(Ql - a2)|| L oo (r . C 3)J \\Zl\\ L 2 (m . y) ||^2|| L 2(9fi;y) • 

Here, as usually, the norm of L°°(T; C 3 ) is 



\w\ 



L°°(T;([ 



£ 

3=1 



W 



(i)| 



lz,°°(r) 



for any vector field w. It is a routine computation to check that, on one hand 

\\K 1 -K 2 \\ Laa{r) <CB(8 c (Cl,Cl)), 
\\V(a 1 -a 2 )\\ Lao(r . >c3) <CB(8c(C^C^), 
\\V(/h-m^<rtc>)<CB(6c(Ch,(%)). 

and on the other hand, 



-1/2 



C^iT) 



+ 



-1/2 -1/2 

-1/2 -1/2 

7i -7 2 



-1/2 

72 



co.^r) 



co.^r) 



c^^r) 



+ 



1/2 



C.^r) 



+ 



1/2 

72 



co.!(r) 



<C 



1/2 1/2 

1/2 1/2 

7i -7 2 



com <CB(8 c (ClCD), 



With all these estimates and Lemma HH in mind, we get 

I ((Qi-Q2)Z 1 \ 3 r 2 ) a \< CB(5 c {ClC 2 v )) ||^|| Bl / 2(8fW INI B v= 



(3ft;y) 



+ 
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+ CB{5 c {Cl,Cl)) (\\Nx<e x \\ THo{t) + \\NxXi\ r \\ TH{r ^ 
x (JWlB^r) + H^i lr H-riY(r) + p2|r|| B i/2 (r ) + IMI Tjf / (r)^ 
hence we deduce the estimate given in the statement. □ 

2.3 Recalling the existence of special solutions 

Let B(0;p) be the open ball centered at the origin O with radius p > and 
such that Q C B{0; p). Sometimes B(0;p) will be denoted by B to simplify 
the notation. Let Eq and po denote the electric and magnetic constants, 
respectively. Extend the coefficients Pj^j defined in Q to functions in E 
-still denoted by jlj, jj~, preserving their smoothness and in such a way that 
Pj — /A),7j — Sq have compact support in B(0;p) (regarding to extensions 
see |24|). Note two simple facts. Firstly, the extensions still satisfy the 
a priori bound and the a priori ellipticity condition in E. Secondly, the 
extensions of the matrices flM]), ( 125]) (125"|) -still denoted by Qj, Q'j, Qj- satisfy 
that u 2 e p Is + Qj, u 2 e p I% + Q'j and u 2 e p I% + Qj have compact support 

in mp^). 

In the following, we state two propositions which were proven in [6] . Their 
proofs are based on ideas from [22], [3], [21] and [T3] . 

Proposition 16 Let 6 be a constant such that —1 < 5 < and let ( G C 3 
be such that ( ■ ( = co 2 e p with 

|C| >C(5,p) J2 \\ uj2£ ^o + qj\\ Lac(B) + Yl \\(u 2 e Voh + Q) k 



B°°(B) ' Z_^i II V~ ""ur-u-o i «W.7|| LO o( S ) 
3=1,2 3,fe=l 



where 



?i = -\*$ -k 2 - \W-DP), q 2 = - l -/\a - k 2 - -^Da-Da). 

Then, there exists a 

Z = e^ x {L + R) 

solution of (-A/ 8 + Q)Z = in E, with Z\ Q E H 2 (Vt; y), 



"a 



1/2 1/2 D 

C • B 
\ uel /2 pl /2 A J 
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for A, B constant complex vector fields, and R satisfying 

C(S,p) 



\ R hy^ id 



\L\J2 ||(w%/i 4 + <5) 



j,k=i 



3 II L°°(B) 



Furthermore, Y = (P — W l )Z is solution for (P + W)Y = in E and it reads 

Y = ( p l l 2 H l | j^E* )* 
with E, H solution for ([1]) in E. 
The norm in the proposition is 



^ 



:i+\x\r\f\ux. 



Again y is meanless, it just stands to remark the form of the elements for 
which the norms are taken. 

Proposition 17 Let ( E C 3 be such that C ' C = w2 ^o/ i o with 



\(\>C(p)J2\\(" 2 ^ i 8 + Q) k j 
j,k=l 



L°°(S) 



Then, there exists a 



Y = e K ' x {M + S) 
solution for the equation (P + W*)Y = in E, with Y\ a E H 1 ^; y), 



M =K\ 



( C-A \ 

-Cxi 

(■B 

V Cxi? ) 



for A, B constant complex vector fields, and S satisfying 

8 



icii < 9ip 2 3. 

Pllia(fl;y)^ |£| 



/m 



L°°{B) 



+ \\w. 



J \\L°°(Q) I ' 
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2.4 Proof of the stability 

The main idea in this final part goes back to pQ. 

Let /ii,7i and /*2,72 be two pairs of coefficients as in Theorem [1] and 
choose 

Ci = ~\i + i (r 2 + !^) V2 m + (r 2 + u/Wo) V2 Tfi, (39) 

C 2 = ^-^(r 2 + ^) 1/ \ 1 + (r 2 + c% / , ) 1/ % 2 , (40) 

with r > 1 a free parameter controlling the size of jCil an d IC2I, where £, 771, 772 
constant vector fields satisfying |r/i| = |r/ 2 | = I, r/i ■ r] 2 = 0, r/j ■ ^ = for 
j = 1,2 and C, j^ e^. More precisely, if £ reads in the coordinates £ as 

£=U (1) e (2) £ (3) )', 

we choose 

with |£'| 2 = (£«) 2 + (£( 2 )) 2 . Observe that Ci - C2 = ~Z and 

£r = 4= + ^= + 0(O, £ = -^ + ^ + ^(r- 1 ). 

Kil \/2 v/2 K2I \/2 y/2 

We now choose other euclidean coordinates J 7 by fixing the following or- 
thonormal basis of M 3 : 

f2 = w\\ e(2) ) ' /a = ( ° ' = e3 ' fl = hxfs 

Here the vectors / x , f 2 , fa are expressed in the coordinates E. In these new 
coordinates £, 771 and 772 read as 

i=\ in 1 , »?i 

Obviously, the metric e in these coordinates is still the identity matrix. 
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Therefore, £1 and ( 2 read in these coordinates T as 



Ci 



, (fZ + k 2 ) 1 ' 2 ^ 



(3) 






<2 



m _ ( r 2 + fc2)l/2S| 

V ? + (- 2 + k rm ) 



where k 2 = uj 2 EqPq. 

Consider Z\ = e % ^ vx (Li + Ri),Y\ the solution stated in Proposition [TBI 
corresponding to the pair /ii,7i with |(j| > C(p,M). Recall that 



Ii 



Id 



/ Ci • A, \ 

1/2 1/2 D 

ue ' p ' B l 
Ci-B 1 

\ 1/2 1/2 „ I 



\R\\ 



L^n-y) 



< 



c{p,n,M) 



Additionally, consider Y 2 = e lC - 2 ' x (M2 + S 2 ) the solutions stated in Proposition 
[T71 corresponding to p 2 ,l2 with |£ 2 | > C(p,M). Also recall that 



M, 



IC2 



/ C2-^2 \ 

-C2 x M 

Ql • B 2 

V C 2 xS 2 j 



I Q II <r 

l°2|| L 2(Q ;: y;) ^ 



C(p,fi,M) 

K2I 



Before plugging Zi and 9^ — Y 2 — Y 2 into the estimate given in Proposition 
[To"j we establish a quantitative version of the Riemann-Lebesgue lemma. 



1 /"!©"> 



Lemma 18 Let r be a positive parameter, q e L 



u q (r) := sup ||g — ^(- — 2/)|lii(R«) 

|j)|<r 



and 



Consider 



GC 



l/"in>n. 



then for any < d < 1 one has 



g m ^g 



r Wy 



^ / n ^ -1 l + |V0(x;r)| .. .. 

<w g (/) + CV SU P i I iT7^„.._M2 H^ll^ 1 ^") 



xm ™ 1 + |V0(x;r 
Proof: Take <p G C °°(lR n ; R+) such that |M| L i (Rn) = 1 with 
suppy? C {x E W 1 : |a;| < 1} 
and denote ip d = d~ n (p(./d). Then one has that 



e m-;r, q 



'a dx 



j4>{-;-r)( n _ 



(q - (f d * q) dx + 



3 «0(.;t) 



(fi * qdx. 
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On one hand, 

Ik - <p< * g|lri(R») < / <p(v)\\q-q(.'-'(v)\\LLCB»)dy. 

On the other hand, since 

(1 + V0(.; r) ■ D)e**<- JT > = (1 + |V0(.; r)| 2 )e^ (-;r) 
one has integrating by parts 



/ 



3 i0(.;r) 



<p d * qdx 



hence 



i<f>{-\r) 



(ft * qdx 



< Cd l sup 



T) f l-V<K.;r)-D 

l + |V0(.;r)|2 

l+|V</>(x;r) 



(^ * g) ete, 



eR ^l + |V0(a;;r)|2 



Mll,i(R™) 



D 



Recall that fij,jj £ H 2+S (£l) with < s < 1/2. In particular, d a ^j,d ar )j 
are in H S (Q) for < |a| < 2. Moreover, the extension by zero allows to 
identify #*(fi) = #*(fi) for -1/2 < t < 1/2 (see [25]). Hence, lnid 01 ^) and 
lQ(9 a 7 J ) are in Hq(Q) which implies 



sup Uln^j - ( 1 n<9 a 7 i )(. - y , Ml/] 

[y|<r 

sup \\lnd a Hj - (lad a ^j)(. - 



L l 



<r 



<Cr s , 
<Cr s . 



(41) 
(42) 



This property can be found in [21 

Now we plugging Z\ and fr 2 = ^2 — Y 2 into the estimate given in Propo- 
sition [TS] with different choices of Aj,Bj. 

Choosing B\ = B 2 = and Ax, A 2 such that 



Vi V2 , , 






one gets, when r becomes large, that 

{{Qi-Q*)z 1 \Y 2 ) Q = 

= J e~ l ^ x ( ^A(ai - a 2 ) + J(Va r Vai - Va 2 -Va 2 ) + {k\ - k\)\ dV 

+ 0((r 2 + |e| 2 )- 1/2 ), 
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where the implicit constant is C(p,Q,M). In addition, by Lemma [TBI with 



4>{x 



and (JUJ), 012]) one has 



r ) = _|£'| x 2 + 2(r 2 + ^2)1/2^3 



(Qi - Q 2 )^i 



F, 



0\d s + 



kl 



■^1 



o v (|£| 2 +|£| 2 |£'l 2 + 40- 2 + * 2 )l£1 2 ) 1/2 

Choosing A\ = A 2 = and JBi, 5 2 such that 
one gets, when r becomes large, that 

{{Qi-Q2)Z 1 \Y 2 ) = 
= Je-*-* Qa(/3 x - fa) + i(VA-VA - V/3 2 -V/3 2 ) + {k\ - k\)\ dV 

+ 0((r 2 + |£| 2 r 1/2 ), 

where the implicit constant is C(p,Q,M). Again, by Lemma ITBl and (|41|) . 
one has 



(Qi - Q 2 )Zi 

Writing 



Yn 



0[d s + 



'-'M 



(le| 2 +le| 2 M 2 + 4(r2 + P)|e|2)i/2y 



/ = lo ( -A(«! - a 2 ) + -(Va r Vai - Va 2 -Va 2 ) + (k 2 - k 2 ) 
= l n ( ^A(/3i - &) + -(V&-V& - V/3 2 -V/3 2 ) + ( K 2 - k 2 ; 



where 1q is the indicator function of Q. By Proposition [15] and the properties 
of the special solutions, there exist three constants c = c(Q), C = C(p, Q, M) 
and C = C'(p, M) such that, for any r > C one has 

1/(01 + \g{i)\ < C (B{5 c {ClCl)y^\ 2 ) 1/2 + (r 2 + |e| 2 )~ 1/a 



+ <T + 



(iei 2 + iei 2 ie? + 4(r2 + p)ie 



2U/2 
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Note that, for R > 1, one has 



+ \\9\\h- he) = [ a + itfr'imoi' + mi 2 )^ 

J\Z\<R 



2 2 

if-!(E) 



+ / (i + ierwuof+imi 2 )^ 

< C {B{5 c {C 1 ,C 2 ))e< R+ ^ + r- 1 + d s ) 2 [ (l + \r\ 2 y l r 2 dr 

Jo 

+C J^J 1 + l ^ yl Kl 2 + lePlfl 2 + 4(r 2 + F)|£'| 2 ^ 
+ (I + R 2 )' 1 (\\f WUn) + \\9\\l H n) 
Lemma 19 One has that 

-2| t |2 ft 



Lj 1 + le ' 2rl Kl 2 + \m? + 4(r 2 + P)|f| 2 di ~ C 



2-' 



<^t 



Proof: Since {(GK 3 : |f | < i?} C {£ 6 M 3 : |f'| < R, |£ (3) | < i?}, the 
integral in the statement is bounded by 

w\< r j\^\<r iei 2 + iei 2 iei 2 + 4(r2 + F)iei^ *■ 

Changing to cylindrical coordinates it is enough to study 

r{r 2 + t 2 ) 

y 1 ^ '" "^ fc ^ 9 I .9 , / 9 
'[0,R]x[0,R] 

One has 



j( Rr ). = [ fi + r 2 + t 2wi r ( r +t J dtdr 



I{R,t)< / (i + r 2 + £ 2 ^ 



r 



2 +t2) l/2 



[0,RM0,R] ( r 2 +t 2 + (r 2 + t 2 )r2) l/2 

r(r 2 + t 2 ) x / 2 
((r 2 + £ 2 )r 2 + 4(r 2 + A; 2 )r 2 ) 1 A 



dtdr 



,, 2 2x i 1 r^ + t 2 ) 1 / 2 

[0jR]xM l + + J (l+r 2 )V2( (r 2 + t 2 )r 2 + 4^2 + ^2)1/2 

<v^ / (l + r' + t 2 )' 1 , - - f — 9 . u „ dtdr 

Jo<r<t<R (t 2 r 2 + 4(r 2 + A; 2 )r 2 )V2 

/" r 2 

+V2 / (l+r^ + t 2 )- 1 — - — -———dtdr 

Jo<t<r<R {r^ + 4(r 2 + /c 2 )r 2 )V2 
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= ^Lr] (t 2 + 4(r2 + F))V2(i +t 2)i/ 2 arctan ((l + t2 ) i/ 2 ) dt < C ~- 

D 
Therefore, 

II/IIh-i(E) + Ibll^CE) < C {B{5 c {Cl,Cl))e< R ^ + t- 1 !? 1 / 2 

+^ s J R 1 / 2 + J R 1 /2^-i r -i/2 + jR -i^ 

Now we choose R in such a way that d s R 1 ^ 2 + i? 1 / 2 ^ _1 r -1 / 2 behaves as 
i? _1 , that is, 

,2/3 1/3 

D _ _ I 

(l + «f 1 + s r 1 /2)2/3' 
hence 



H-!(E) 



blln-iCE) < C fs(5 c (C^,C 2 ))e CT + (V + ^^ J 



Choosing r — d 2 ( 1+s ) the estimate becomes 

II/L- 1( e) + IMI*-i(B) < c (s(fc(c^,^))e^- 2(1+i 

On the other hand, the a priori bound was chosen to have 

Ih-(O) + IMI/r«(n) ^ C '( M )' 



for < s < 1/2. Finally, interpolation theory ensures the existence of two 
constants C = C'(p,M) and C = C(p,Q,M,uj) such that, for any d < C , 
the following estimate holds 



L 2 (fl) 



+ IMI^n) < c{B{5 c {ClCl))e^- 2{1+s) + £%)° (43) 



with = -e + (l-0)s. 

The idea now is to transfer this estimate from /, g to the difference of 
the coefficients jX\ — \X2 and 71 — 72. This can be accomplished by using the 
following Carleman estimate. 

There exists a positive constant C(Q) such that, for all h < 1 and any 
function smooth enough, the following estimate holds 

/ille v/ ^H 2 + /i 3 lle^V^II 2 < 

'Mr ^IIl 2 ^) ^ "" Ir v< ^llL 2 (n ; c 3 ) - 

< C (V He^A^II 2 + /ille^e!>l| 2 + /i 3 lle^V^II 2 ^ 
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where ip = \j1\x — Xo| 2 with x ^ Q. The constant here depends on the 
distance from x to Q and on the diameter of Q. A Carleman estimate of 
this type can be found in [9]. 
A simple computation give: 



f 1 ~-l/2 

/ = ln7i 



A/-1/2 -l/2x . /-1/2 ~l/2s . ,-1/2 -1/2 

A(7i -7 2 ) + 9/(7i / -7 2 )+Wi - " 



1 —1/2 

9 = lnA*i 



•~l/2 -1/2 



-1/2 ~l/2\ 



^2 

-1/2 ~l/2x, 



A f ~ i. z ~i-/-\ i t~t z ~J-/^\ i t~- - ~ J- /^\ 

A(/V ~^2 ) + ?s(/V -^2 )+JWY -7 2 ) 



where 
9/ = " 



. 2-1/2,-1/2- . -1/2- x \ 2- -1/2,-1/2 . ~l/2x 

+ w 7l (7l Ml + 72 M2) , Pf = -U 7l7 2 Ol +/*2 )) 



A ~V 2 

-1/2 
72 

A/2 2 2-1/2/-1/2- -1/2- x 

— TFr + W Ml Ol 71 + M 2 72) 

M 2 



2,^,-1/2^1/2 ,- 1/2, 



2- ~i.1zf~i.1z . ~i/z\ 
p g = -U MlM 2 (7l +72 )• 



Note that, thanks to the a priori bound, we have the following differential 
inequalities: 

| A /~l/2 ~l/2xi ^ nC\/t\(\ f\ 1 I- 1 / 2 - 1 / 2 ! 1 l~l/2 ~l/2ix 

l A (7i -7 2 )l < C(M)(|/| + |7i' -j 2 ' I + |Mi -fh \h 

I a/ -1/2 -1/2m ^ nf\/r\(\ I 1 I- 1 / 2 - 1 / 2 ! , l~l/2 ~l/2ix 

|A(/V -// 2 )| <C(M)(|y| + |7i -72 I + ]Mi -IH. I)- 

~ 1/2 1/2 

In order to simplify the notation, we shall write <pi = 7/ — 7 2 and 



1/2 -1/2 

Pi 



^2 : : f~i - /V • By the differential inequalities written above and the 
Carleman estimate, one has 



E(* 



^A, 



^IIl 2 (s1) 



+ /i 3 ||e v//l 



v ^L 2(n . c3) ) < 



3=-,- 



< nil V^ f.4 11 <p/h ± ii 2 j_ . ILvA^.H 2 _i_ 7,3 



3 v/^ 



V0, 



■?IIl 2 (<9S1;C 3 ) 



i=i, 2 



-\-u n \ ne j r2 , x + e </ 



|2 
Il 2 (s1) 



L 2 (S1) 



where the constant is C" = C"(Q, M) and (f(x) = l/2|x — Xo| 2 with x ^ fi. 
The terms h 4 || e ^ h 0i|| L 2 m x) with j = 1,2, can be absorbed by the left hand 
side of the inequality. Hence, if d\ = inf{|x — xq\ 2 '■ x G Q} and d 2 = 
sup{|x — x \ 2 : x G Q} we get, for any h < C"(Q, M)" 1 / 3 , that 



h 4 



e dl/h E (* H^ll^(n) + h * HW,-|£ 2( n ;C 3)) < CV*" /h x 
i=i,2 

L 2 (S1) + IML^J + / j ( J 1 W^jW^idQ) + ^ IIWillL2(aO;C3) 



i=i,2 
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But now we can easily estimate 

WMweti) + HV0i|| i2( an;C3) < CB(5 c (C l r ,C 2 )), 

WfaWvyxi) + llV0 2 || L2( an;C3) < CB(6 c {Cl,C*)), 

Il7i - 72|| L 2 (n) + ||V(7i - 72)|lxa ( n ; cs) ^ C (ll^ll^n) + IIWiIIl^cs)) > 

lljSi - /i 2 || L 2 (n) + ||V(/ii - Ab)||, L 2 (n . C 3 ) < C ^||0 2 || L 2 (n) + ||V0 2 || L 2 (n . c3) J . 

The constants above depend on the a priori bounds M. With these inequal- 
ities and estimate (1151) . we obtain 

rf2~ d l 



1 7i - 72 || H i(n) + llAi - ^11 H i (n) ^ ^ 2h B(5 C (C^C£)) 
+Ce^(B(5c(C 1 r ,C*)y<- 2il+e) + ^ 



2s \ i+s 



where d 2 > d h < s < 1/2, C = C{p,Q,M), i < C'{p,M), c = c(fi) and 
h < C"(Q, M) _1//3 . To end up with the estimate given in the statement, it is 
enough to choose the parameter d as 

J- 2 ^ = -ho g B(5 c (Cl,C 2 r )), 

and to note that 

s 2 s 2 

0< 3i7TTr < r 

3 The domain U is partially spherical 

Along this section we assume U to be a suitable partially spherical domain 
and we follow the notation in Definition [10] and Definition [TTJ Furthermore, 
n and v will denote the outward unit normal forms of U and Q, respectively. 

The basic idea in this section is to use the Kelvin transform /C to gen- 
eralize our result on partially flat domain to the case of partially spherical 
domain. To achieve this, we study the behavior of Maxwell's equations and 
the distance 5c under K. 

Note that /C = /C _1 and /C is a conformal transformation from (Q, e) onto 
(17, e): 

4 

JC*e = Are, 



where fC* denotes the pull-back of K. 

Let E = K*E, H = }C*H, jj, = JCfi, and 7 = /O7. The following is the 
transformation law for Maxwell's equations under the Kelvin transform. 
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Lemma 20 One has E,H G H(U; curl) is solution of 

dH + ?W7 *E = dE — ioo/j, *H = 
in U, if and only if, E,H G if (fi; curl) is a solution of 

2 2 

d/T + iw7j^2*.E = dE - «w/x-pJ2 *if = 

in f2. 

Proof: The proof follows easily from 

dJC* V = JC*d v , K*(* v ) = * K . e K% * C eV = c 3/2 ~ k *v- 

Here f] is fc-form and c is an arbitrary positive smooth function. □ 

Lemma 21 Given Uj G H(U; curl) and Vj G H(Q; curl) with j = 1, 2, let us 
consider Uj = JC*Uj G H(Q; curl) and i>j = K,*Vj G if (£/; curl). 

(a) For any z G B 1 / 2 (dU;A 1 TE) one has (*(nAWj)|«) = (*(v A Uj)\w) , 
where w = /C*vjarj G £ 1/2 (<9fi; A*TE) with u G if^C/; A X TE) such that 
v\du — z - Furthermore, 

\\ w \\B 1 /' 2 (dn;A 1 TV) — C lkllBi/2(at/ ; AiTE) • (44) 

For any w G 5 1/2 (9Q; A 1 ! 1 !) one has (*(i/ A Vj)\w) = (*(nAVj)\z), 
where z = JC*u\ du G B 1 ' 2 {dU; A 1 TE) with u G if : (Q; A : TE) such that 
u\do. = w. Moreover, 

W Z \\B 1 / ,2 (dU;A 1 TE) — ^ \\ W \\ BV 2 (dtl-^TK) ' (^) 

(b) For any /i G B 1 / 2 (dU) one has (Div*(n Aiij)|/i) = (Div *(i/ A iij)\g) 
where # = /C*/U G E 1 / 2 ^) with / G H l {U) such that /| a[/ = h. 
Moreover, 

lbll B i/2 (9n ) <C II^Li/a^) • ( 46 ) 

For any g G i? 1 / 2 ((9f2) one has (Div *(v A Vj)\g) = (Div *(n A Vj)\h) 
where h = K*f\ du E B l / 2 (dU) with / G if 1 (ft) such that f\ da = g. 
Moreover, 

ll^llsi/2(a[/) < C 1 1 51151/2 ^n) • (47) 
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(c) The following estimates hold 

||*(nA Ui) - *{n /\u 2 )\\ TH{du) < C \\*(v A Ui) - *{v Au 2 ) 



and 



|*(i/Avi) -*(^Aw 2 )|| TH(an) < C"||*(n Avi) -*(n Ai> 2 ) 



iT_ff(en) 



\TH(dU) 



Proof: The proof of the identities is an immediate consequence of the 
identities stated in the proof of Lemma [201 and the weak definitions of tangen- 
tial trace and surface divergence. Proving (j46p . (1471) is an easy computation 
and (j44"|) . (J4"5l) follow easily in coordinates from (J46]), (J4"T|) and ffTTj) . Finally, 
the estimates in (c) are a consequence of (a), (b), (J3J) and (J3J). D 

Proposition 22 One has that 
where 



cf = c(^,7i;r), c£ = c 



In ' 1 - 



r^jj 



7i;r 



with j — 1,2. 



Proof: Considering Ej,Hj and Ej,Hj as u^ and Vj in the statement of 
Lemma EH this proposition is a consequence of Lemma [20] and the item (c) 
in Lemma I2~T1 □ 

In order to end up with the proof in the case that U is partially spherical, 
it is enough to use Proposition [221 and recall that 



1/^1 - A t 2||^i(( 7 ) < C 

I|7i-72|Ihi(u) <C 



\ifi\-fa) 



-1^(7-72) 
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